#ifdef __cplusplus
extern "C" {
#endif
#include <stdio.h>
#include <string.h>
#include <unistd.h>
#include <math.h>  

#include "strsubs.h" 

#define	CHI_EPSILON     0.000001    /* accuracy of critchi approximation */
#define	CHI_MAX     99999.0         /* maximum chi square value */

#define	LOG_SQRT_PI     0.5723649429247000870717135 /* log (sqrt (pi)) */
#define	I_SQRT_PI       0.5641895835477562869480795 /* 1 / sqrt (pi) */
#define	BIGX           20.0         /* max value to represent exp (x) */
#define	ex(x)             (((x) < -BIGX) ? 0.0 : exp (x))
#define	SQRT_PI       (1.0/I_SQRT_PI)  /* sqrt (pi) */

double medchi(int *cls, int len, int *n0, int *n1, double *kstail)  ;
double ks2(int *cls, int len, int *n0, int *n1, double *kstail)  ;
double probks(double lam) ;

double nordis(double z) ;
double ndens(double val, double mean, double sig) ;
double ntail(double z) ;
double zprob(double p) ;
void setzptable() ;
double z2x2(double *a) ;
double conchi(double *a, int m, int n)  ;
double conchiv(double *a, int m, int n)  ;
double chitest(double *a, double *p, int n) ;

double xlgamma(double x) ;
double psi(double x) ;
double tau(double x) ;
void bernload() ;
double bernum(int x) ;

void mleg(double a1, double a2, double *p, double *lam) ;

double dilog(double x) ;
double li2(double x) ;

double hwstat(double *x) ;

double gammprob(double x, double p, double lam)  ;
double bprob(double p, double a, double b)  ;
double lbeta(double a, double b)  ;
double dirmult(double *pp, int *aa, int len) ;
double dawson(double t) ;

double binomtail(int n, int t, double p, char c) ;
double binlogtail(int n, int t, double p, char c) ;
void genbin(double *a, int n, double p) ;
void genlogbin(double *a, int n, double p) ;
int ifirstgt(int val, int *tab, int n)  ;
int firstgt(double val, double *tab, int n)  ;

void cinterp(double val, double x0, double x1, 
  double f0, double f0p, double f1, double f1p, double *fv, double *fvp)  ;
int firstgtx(double val, double *tab, int n) ;

double rtlchsq(int df, double z) ;
double critchi(int df, double z) ;
double rtlf(int df1, int df2, double f) ;

double ltlg(double a, double x)  ;
double rtlg(double a, double x)  ;

double twdens(double twstat)  ;
double twtail(double twstat)  ;
double twtailx(double twstat)  ;
double twdensx(double twstat)  ;
double twnorm(double lam, double p, double n)  ; 
void twfree() ;
int settwxtable(char *table) ;
void gettw(double x,  double *tailp, double *densp)   ;
double dotwcalc(double *lambda, int m, double *ptw, double *pzn, double *pzvar, int minm) ;
int numgtz(double *a, int n) ;

double betaix(double a, double b, double lo, double hi)  ;
double betai(double a, double b, double x) ;
void bpars(double *a, double *b, double mean, double var)  ;
void bmoments(double a, double b, double *mean, double *var)  ;
double unbiasedest(int *ndx, int ndsize, int **counts) ;
void weightjack(double *est, double *sig, double mean, double *jmean, double *jwt, int g)  ;
int modehprob(int n, int a, int m) ;
void calcfc(double *c, int n, double rho) ;
void circconv(double *xout, double *xa, double *xb, int n) ;

double bino(int a, int b) ;
void setbino(int maxbco) ;
void destroy_bino() ; 
#ifdef __cplusplus
}
#endif
